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NOTATION 


plate length in x-direction 

plate length in y-direction 

plate thickness 

lateral pressure 

Vertical displacement of points of the middle surface 
membrane stress parallel to x, and y axes 

membrane shearing stress 


extreme-fiber bending and shearing stress parallel to 
x and y axes 


bending shearing stress 


bending moments per unit length of sections of a plate 
perpendicular to x and y axes, respectively 


unit elongations 1n x and y direction 

shearing strain component 

modulus of elasticity in tension and compression 
poisson s ratio 

flexure rigidity of the plate 

stress-funciion 


Subscripts m.n.m,n.r, Ss. p gq, k, tand § represent 
integers. 
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ee menOduction 

Plates, generally speaking, deflect as nondevelopable surface. 
As a consequence, the increase of deflections beyond a certain magnitude 
is accompanied by the so-called membrane stresses acting in the plane of 
the plate. In order to describe the behavior of the plate at deflections 
equal to or exceeding its thickness but still small compared to its lateral 
dimensions, Von Karman derived the large-deflection theory of isotropic 
plate by assuming that the external loads are resisted by both bending ac- 
tion and forces in the middle plane of the plate. 

Von Karman'‘s fundamental equations for large-deflection theory 
has been solved by many authors, but most of them are of an approximate 
nature. The only theoretically exact solutions were given by S. Levy. But 
he gave the results for square plate only since his theoretically exact so- 
lutions are extremely cumbersome. So, up to the time being, we still do 
not know exactly the behaviors of a rectangular plate with high aspect 
ratio at large deflection. As it is well known in designs of thin plates 
that bend under lateral loading, formulas based on small deflection theory 
which neglects stretching and shearing in the middle surface may still be 
used with sufficient accuracy for deflections up to fifty per cent of the 
plate thickness. This is true for a square plate or a rectangular plate with 
low aSpect ratio. If a rectangular plate with high aspect ratio is being 
considered, we cannot make sure whether it is true or not, since no exact 
results have been given for such a rectangular plate at large deflection. 


It is the purpose of this paper to get an insight into the behaviors of 





plates with high aspect ratio by means of S. Levy's theoretical exact 
method. 
It should be mentioned here that this paper has been restricted 
eiauto, 
a) all edges simply supported 
b) lateral uniform pressure, total edge load zero 
c) aspect ratio es ao 
This shows that the areas for further research have been by no means ex- 
hausted. 
The choosing of aspect ratio 6 = 5 is not arbitrary, since this 
value is the upper limitation in the ship structure. Moreover, beyond this 
value, the behavior of the plate is nearly as a strip. And itis, of course, 


out of the purpose of this paper if the value of the aspect ratio is small. 


2. Fundamental equations for the deformation of thin plates 

Consider a rectangular isotropic elastic plate having the thick- 
ness h andrigidity D. Let the origin of a cartesian coordinate system 
x, y, z be located at a corner of the middle plane of the plate, the axes 
x and y being directed along the edges of the plate of length a andb, 
respectively. The plate is acted on by an arbitrary lateral load p(x, y), 
which provokes a displacement field in the plate characterized by the 
components u, v, w, the last one representing the deflection of the 
plate. The fundamental equations governing the deflection of thin plates 
were developed by Von Karman. They were given by Timoshenko [3]* in 
essentially the following forms: 


*Numbers in the brackets, []|, refer to references at the end of 
this paper. 
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where the membrane stresses are 
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and the extreme-fiber bending and shearing stresses are 
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3. General solution for simply supported rectangular plate 
In case of simply supported rectangular plate, the boundary 


conditions referred to deflection and moment are 


Oe TX, — D538 t+ $72) =e hin HHO, O (6a) 
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One of the boundary conditions referred to the stress function F 


(6b) 


is that there are no shearing stresses along the edges. This boundary 
condition can be written as 


KH = O at the boundaries; (7) 
another is that the boundaries remain straight but may displace perpen- 
dicular to their lengths. [t should be noted that, if this is done, the 
distribution of the in-plane stresses perpendicular to the boundaries 
cannot be specified, and we can merely state that the integral of the in- 
plane stress along the boundary must be equal to the total applied load 
along the boundary. 


An exact solution of this problem can be established by starting 


from the simultaneous equations (1) and (2). 
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The deflection of the plate may be taken in the Navier form 


i) ’ mix \ se, 
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the boundary conditions (6) thus being satisfied by any, yet unknown, 
values of the coefficients w,, ,. The given lateral load may be expanded 


in a double Fourier series 
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are constants denoting the total tension load applied on the sides x =o, 
a and y=o, b, respectively. Substituting the expressions (8) and (10) 
into (1) we get the following relation between the coefficients of both 


series: 
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p,q’ are to be determined by the val- 


ues of the coefficients We The sum includes all products for which 
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Inserting the expressions (8), (9) and (10) into equation (2), we 


can establish the relation between the deflections, the stress functions and 


the lateral pressure. The general equation given by midi is as follows: 
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4. Specific solution for rectangular plate with high aspect ratio under 
uniform lateral pressure, total edge load zero 

a. General: In this special case, the total edge load in the x- 
direction kb and in the y-direction p, is zero, and the uniform lateral 
pressure 1S constant. The expansion of this pressure P in a Fourier ser- 
ies as shown in equation (9) gives pressure coefficients ae = ot (4Y p. 
In the solution of a rectangular plate with such a high aspect ratio as = on 
we assume that, in the deflection coefficients War m= 1 and n=15 
are accuracy enough. So we consider w) ),W) 3,W] 5,W] 7,W] 9,Wj 11> 
W1,13 and w) 75 only. All other terms assumed to be equal to zero. 


Under this assumption, a general equation for the coefficients in the 


stress-function can be obtained from equation (12) as follows: 
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Equation (13) represents a doubly infinite family of equations. 


In this special case, equation (13) can be reduced as 
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where S=1, 3, 5, ---- 15, representing the first eight simultaneous 
equations of the family. 
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In each of the equations of the family the coefficients for the 
stress function, tog , may be replaced by their values as obtained from 
equation (14). The resulting equations will involve the lateral pressure 
p and the cubes of deflection coefficients Ween . oolving these eight 
simultaneous nonlinear equations, we obtain the values of the eight de- 
flection coefficients, Wi] a] oe MG The results are given in 
Table (1) and shown in Fig. (1, 2). The straight lines give the correspond- 
ing solution by the linear theory [3]. 


The equation of the deflection profile for the center line parallel 


to the long sides of the plate can be written as: 
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a part of the results are given in Table (2), and three curves are shown in 


Pg (3). 


b. Membrane stresses: From equations (3) and (10) we can 


derive the equations for the membrane stresses as follows: 
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Substituting the values of deflection coefficients w given 


m,n 


in Table (1) into equation (14), and then substituting equation (14) into 
equation (17), we obtain all the results, a part of which are given in 
Tables(3-7). Three of each group are shown in Figs. (4, 5), and the mem- 


brane stress curves for some specified points are shown in Fig. (8). 


c. Bending stresses: From equations (5) and (8) we can derive 


the equations for the bending stresses as follows: 
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Substituting the values of deflection coefficients Wie given in Table (1) 
into equation (18), we obtain all the results, a part of which are given in 
Tables (8-11). Three of each group are shown in Figs. (6-7). Bending stress 


curves for some specified point and maximum bending stresses curves are 


shown in Fig. (9). 


d. Stress intensity: For determining failure of ductile materials 
having a pronounced yielding point stress Oy p in simple tension, octa- 


hedral failure criterion seems to be the most significant and applicable 


theory [2]. 
Define “y= te + fe (19a) 
oy = 6, +0 (19b) 
y, +/ 
Cyy= Ty t Oy (19c) 


Octahedral failure criterion can be written as 


oop = [c+ é, = by by + 3b (20) 


From the results we have obtained, we can see that the values 
of the maximum bending stresses are much more larger than that of the 
maximum membrane stresses, and the maximum bending stresses occur on 
the center line parallel to the long sides of the plate where the shearing 
stresses are zero. By introducing the dimensionless notation, equation 


(20) becomes 
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The results are shown in Fig. (10). 


9. Conclusion and discussion 

a. From curve (2) we can see that, for design purpose, even in 
the case of rectangular plate with high aspect ratio, the linear theory may 
still be used provided that the deflections at center are small compared 
with the thickness. But it should be noted that the error in deflection is 
about 15 per cent when the deflection at center is equal to 50 per cent of 
the plate thickness. Comparing with the results of square plate, we can 
see that the error is increased as the aspect ratio of the plate is increased. 

b. It is interesting to point out here that, as shown in Fig. (3), 
for a rectangular plate with high aspect ratio, the maximum deflection 
shows up at the center of the plate when the deflection is small but not at 
the center when the deflection is large. This is probably due to the high 
compressive membrane stresses Oy occurring near the corners of the 
plate. It is evident that this compressive membrane stress is the result 
of the boundary condition that the edges remain straight but may displace 
perpendicular to their length. 

c. As we can see from Fig. (8) that the membrane stresses are 
negligible provided that the deflections are small compared with the plate 
thickness, and become more and more significant as the deflections of the 
plate increase, and that the membrane stresses in the x-direction are al- 
ways larger than that in the y-direction. It is now clear that the error in 
deflection generated by using small deflection theory, up to 15 per cent 


for deflection equal to 50 per cent of plate thickness, is mainly due to the 
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neglecting of membrane action in x-direction. 

d. The fact that the absolute value of the membrane stresses in 
y-direction at the corners of the plate are just equal to that at the center 
of the short sides of the plate is, of course, due to the one term approxi- 
mation in that direction. 

e. From Fig. (9) we can see that the maximum bending stresses 
in the x-direction are prominent when the lateral load, as well as deflec- 
tion is small, and the maximum bending stresses in the y-direction be- 
come prominent when the lateral load is large. This is due to the deflec- 
tion hump occurring as the lateral load is large enough. 

f. The distribution of all the stresses along the y-direction 
changes very rapidly near the ends of the plate and is nearly equal to 
constant in the middle part, about 60 per cent of the plate length. This 
indicates that the middle part of the plate behaves like the case of 
cylindrical bending. 

g. Fig. (10) shows that the lateral load at which the yield 
stresses occur is dependent upon the yielding point stress of the material 
and the width-thickness ratio (F) . Once we know the load, we can find 
the deflection from Fig. (2). Example: 
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Table 2. Values of deflection profile, Fz . along x= 
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